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I. INTRODUCTION 

Inclusive lepton scattering has for many decades been 
the most important tool with which to probe the inter¬ 
nal quark and gluon (or parton) structure of nucleons and 
nuclei. Structure functions extracted from inclusive deep- 
inelastic scattering (DIS) experiments display the central 
features of quantum chromodynamics (QCD) — asymp¬ 
totic freedom at short distances (via structure function 
scaling and its violation) and confinement at large dis¬ 
tance scales (via parton momentum distributions). 

Since the late 1960s, DIS experiments have yielded an 
impressive data set that maps nucleon structure functions 
over several orders of magnitude in the Bjorken scaling 
variable, x, and the squared four-momentum transfer, 
Q^. These data, supplemented by cross sections from 
hadronic collisions and other high-energy processes, have 
enabled a detailed picture of the parton distribution func¬ 
tions (PDFs) of the nucleon which are extracted through 
global QCD analyses (see Ref. [T] and references therein). 

At lower energies, where nonperturbative quark-gluon 
interactions are important and the inclusive lepton- 
nucleon cross section is dominated by nucleon resonances, 
the structure functions reveal another intriguing feature 
of QCD, namely, quark-hadron duality. Here, the low 
energy cross section, when averaged over appropriate en¬ 
ergy intervals, is found to resemble the high energy result, 
whose dependence is described by perturbative QCD. 
In this context, quark-hadron duality provides a unique 
perspective on the relationship between confinement and 
asymptotic freedom, and establishes a critical link be¬ 
tween the perturbative and nonperturbative regimes of 
QCD. 

In the framework of QCD, quark-hadron duality can 
be formally interpreted in terms of structure function mo¬ 
ments [5]. From the operator product expansion (OPE), 
the moments can be expressed as a series in l/Q^, with 
coefficients given by matrix elements of local quark-gluon 
operators of a given twist. The leading (twist 2) term 
corresponds to scattering from a single parton, while 
higher twist terms correspond to involving multi-quark 
and quark-gluon interactions. Since at low the res¬ 
onance region makes a significant contribution to the 
structure function moments, one might expect a strong 
dependence of the low-Q^ moments arising from the 
higher twist terms of the OPE. In practice, however, the 
similarity of the structure function moments at low 
with the moments extracted at high energies suggests the 


dominance of the leading twist contribution, with the 
higher twist, multi-parton contributions playing a rela¬ 
tively minor role or that the large higher twists cancel 
each other. 

This non-trivial relationship between the low-energy 
cross section and its deep-inelastic counterpart was first 
observed by Bloom and Gilman [5] in the early DIS mea¬ 
surements that were instrumental in establishing struc¬ 
ture function scaling. More recently, the availability of 
extensive, precise structure function data from Jefferson 
Lab, over a wide range of kinematics, has opened up the 
possibility for in-depth studies of quark-hadron duality. 
Duality has now been observed in the proton F 2 and 
structure functions gHZl, the nuclear structure function 
F 2 [H], the spin-dependent gi structure functions of the 
proton and ^He [9l[T0], the individual helicity-1/2 and 
3/2 virtual photoproduction cross sections for the pro¬ 
ton [TT] . 

To establish the dynamical origin of quark-hadron du¬ 
ality one must also study the neutron. Four-quark higher 
twists contributions suggest that duality in the proton 
could arise from accidental cancellations between quark 
charges, which would not occur for the neutron |12j . Un¬ 
fortunately, the absence of high-density free neutron tar¬ 
gets means that essentially all information on the struc¬ 
ture functions of the neutron has had to be derived from 
measurements on deuterium. Typically, the deuterium 
data are corrected for Fermi smearing and other nuclear 
effects HMZl, which introduces an element of model de¬ 
pendence into the procedure. This is particularly prob¬ 
lematic in the nucleon resonance region, where Fermi mo¬ 
tion leads to significant smearing of the resonant struc¬ 
tures. The existence of duality in the neutron F 2 struc¬ 
ture function was suggested recently m in an analy¬ 
sis which used an iterative deconvolution method |18j to 
extract neutron resonance spectra from inclusive proton 
and deuteron F 2 data . A direct confirmation of duality 
in the neutron, however, was not possible. 

Recently, a new experimental technique, based on 
spectator nucleon tagging m, has been used to extract 
the free neutron F 2 structure function m- By detecting 
low-momentum protons at backward angles in electron 
deuteron scattering, the BONuS experiment at Jeffer¬ 
son Lab measured F^ in both the resonance and DIS 
regions, with minimal uncertainty from nuclear smear¬ 
ing and rescattering corrections [21]. In the present pa¬ 
per, we use the BONuS data to quantify for the first 
time the degree to which duality holds for the F 2 struc- 
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ture function of the free neutron. Because the results 
reported here use data from an experimentally-isolated 
neutron target, one can expect significantly reduced sys¬ 
tematic uncertainties compared with those in the model- 
dependent analysis of inclusive deuterium data m- 

For the theoretical analysis of duality we use the 
method of truncated structure function moments devel¬ 
oped by Psaker et al. [52]. Here, the n-th truncated mo¬ 
ment of the F 2 structure function is defined as 

Mn{Xrnin,Xma.^,Q'^) = dx (x, , (1) 


where the integration over x is restricted to an interval 
between Xmin and Xmax- This method avoids extrapola¬ 
tion of the integrand into poorly mapped kinematic re¬ 
gions, and is particularly suited for the study of duality 
where an x region can be defined by a resonance width 
around an invariant mass = M^ + Q^{l — x)/x, where 
M is the nucleon mass. As the position of the resonance 
peak varies with x for different values, the values for 
Xniin and Xniax evolve to correspond to the appropriate 
invariant mass squared region. For the BONuS data, four 
ranges in W'^ were considered, corresponding to the three 
prominent resonance regions (1.3 < < 1.9 GeV^ for 

the first or A resonance region, 1.9 < IF^ < 2.5 GeV^ for 
the second resonance region, and 2.5 < IF^ <3.1 GeV^ 
for the third resonance region), as well as the combined 
resonance region (1.3 < IF^ < 4 GeV^). Results on the 
lowest three non-trivial moments are reported and com¬ 
pared with recent global PDF parametrizations, as well 
as with previous model-dependent data analyses. 

In Sec. |ll| we review the BONuS experiment and the 
results for the neutron F 2 structure function. The anal¬ 
ysis of the truncated moments is discussed in Sec. m 
together with the implications for duality and its viola¬ 
tion. Finally, in Sec. IV we summarize our results and 
discuss their wider implications. 


II. THE BONUS EXPERIMENT 

The results reported here rely on a novel experimen¬ 
tal technique aimed at eliminating or substantially re¬ 
ducing the theoretical uncertainties involved in extract¬ 
ing neutron data from nuclear targets. The BONuS 
(Barely Off-shell Nucleon Structure) experiment [T5H5T| 
used a Radial Time Projection Ghamber (RTPC) to de¬ 
tect backward-moving, low momentum spectator protons 
produced in electron-deuterium scattering in conjunction 
with electrons detected using GLAS [55] in Hall B at Jef¬ 
ferson Lab. By tagging low momentum backward moving 
spectator protons one minimizes final state interactions 
[551 - 155] and ensures that the neutron is just barely off- 
shell [5D] . Additionally, Fermi smearing effects are essen¬ 
tially eliminated. 

The BONuS experiment ran in 2005 and acquired 
electron-deuteron scattering data at two electron beam 
energies: 4.223, and 5.262 GeV. The RTPG consists of 



FIG. 1. Kinematic coverage of the BONuS data. The lines 
represent the IF^ thresholds for the four resonance mass re¬ 
gions mentioned in the text. 


three layers of gas electron multipliers surrounding a thin, 
pressurized gas deuterium target which can detect spec¬ 
tator protons with momenta as low as 70 MeV/c. The 
scattered electron was detected by CL AS. The experi¬ 
ment and data analysis are described in detail in Ell- 
Ratios of neutron to proton F 2 structure functions and 
the neutron F 2 structure function by itself were extracted 
over a wide kinematic range and for proton spectator mo¬ 
menta between 70 and 100 MeV/c. The total systematic 
uncertainty in the neutron structure function extracted is 
8.7% [5T]. Additionally, there is an overall 10% scale un¬ 
certainty due to cross normalization of the BONuS data 
to existing parameterizations. 

The kinematic coverage shown in Figure]^ (4.223 and 
5.262 GeV results combined) extends from the quasielas¬ 
tic peak to the deep inelastic region corresponding to 
final state invariant masses of IF^ > 5 GeV^. The curves 
shown represent the IF^ thresholds for the resonance 
mass regions mentioned above. Typical F 2 results for 
= 1.2 and 2.4 GeV^ are shown in Figure 

The open/closed symbols correspond to data from 4.2 
and 5.2 GeV electron beam energies, respectively. Pre¬ 
dictions, with and without higher twist effects, of the 
ABKM QGD fit [57] are also shown. Qualitatively, one 
can see evidence here for quark-hadron duality in that 
the curves generally go through the resonance data. 

The data used in this analysis had spectator angles 
with respect to the momentum transfer greater than 100° 
and momenta between 70 and 100 MeV/c. 


III. TRUNCATED MOMENTS AND LOCAL 
QUARK-HADRON DUALITY 

Since Q^, x, and W'^ are not independent of each other, 
a range in IF^ at fixed implies a corresponding range 
in X. This allows for a straightforward integration of 




3 



30 

20 

10 


10 


first 

• BoNuS 
* ^ABKM 

1 second 

□ 

• i 

• 

• 

• 

• 

• 

• ,• 

• 

B 

• 

• 

• 

■ • • . 

third 

total 

t 

t 

ff 

• 

• 

• 

% 

a 

• 

a 

1 

• 

• 


30 

20 

10 


1 Z 3 4 1 2 , 3,4 

[GeV^] [GeV"] 


80 

60 

40 

20 


FIG. 2. Typical neutron structure function data from the 
BONuS experiment at = 1.2 GeV^ (top panel) and = 
2.4 GeV^ (bottom panel). The open circles represent data 
obtained with a beam energy of = 4.2 GeV, while the 
closed circles were obtained with E = 5.2 GeV beam energy. 
The curves shown are the ABKM DIS parameterizations m\ 
with (solid line) and without (dashed line) higher twist effects 
and target mass corrections. 


FIG. 3. The second neutron truncated moments M 2 , as a 
function of Q^. The closed circles represent the moments 
obtained from the BONuS data, while the blue rectangles are 
the moments obtained from the ABKM parameterization m 
which includes target mass and higher twist corrections. 



1000*M2 

(GeV^) 

first 

second 

third 

whole 

1.00 

31.5 ± 1.1 

16.4 ± 0.4 

14.1 ± 0.3 

76.7 

± 

1.2 

1.20 

23.5 ± 0.5 

15.3 ± 0.3 

13.0 ± 0.3 

67.4 

± 

0.6 

1.40 

17.7 ± 0.4 

13.5 ± 0.2 

12.1 ± 0.3 

57.7 

± 

0.5 

1.70 

12.3 ± 0.3 

10.7 ± 0.2 

10.5 ± 0.2 

46.7 

± 

0.5 

2.00 

8.4 ± 0.2 

8.5 ± 0.2 

9.0 ± 0.2 

38.4 

± 

0.4 

2.40 

5.8 ± 0.2 

6.1 ± 0.1 

7.1 ± 0.3 

29.5 

± 

0.4 

2.90 

3.4 ± O.l 

4.5 ± 0.1 

5.3 ± 0.3 

21.5 

± 

0.4 

3.40 

2.1 ± 0.1 

3.1 ± 0.1 

4.0 ± 0.2 

15.8 

± 

0.3 

4.10 

1.3 ± 0.1 

1.7 ± 0.1 

2.8 ± 0.1 

N/A 




TABLE I. Second order truncated moments of the neutron 
F 2 structure function based on the BONuS data for the three 
resonance regions studied, as well as for the whole region 1.3 < 
< 4 GeV^. For = 4.10 GeV^ the moment for the whole 
region was not computed due to the lack of data in the highest 
region. 


the experimental F 2 structure function data to obtain 
truncated moments. In order to minimize model depen¬ 
dence, the integrals were evaluated based solely on the 
experimentally measured points without using any inter- 
or extrapolating function. The second (n = 2) truncated 
moments, M 2 , obtained from these data are shown in 
Fig. I as a function of , and are listed in Table |l] The 
uncertainties quoted take into account the experimental 
statistical and systematic uncertainties but do not the 
include the 10% scale uncertainty due to cross normal¬ 
ization of the BONuS data. The corresponding higher 
order truncated moments (n = 4 and n = 6) are pro¬ 
vided in Tables |n] and m respectively. 

In order to study local duality, we formed the ratio 
of the truncated moments of F 2 in the resonance region 



1000*M4 

Q" (GeV") 

first 

second 

third 

whole 

1.00 

11.58 ± 0.43 

3.09 

± 

0.08 

1.69 

± 

0.04 

17.49 

± 

0.44 

1.20 

9.80 

± 

0.21 

3.51 

± 

0.06 

1.95 

± 

0.04 

16.78 

± 

0.22 

1.40 

8.11 

± 

0.17 

3.60 

± 

0.06 

2.17 

± 

0.04 

15.61 

± 

0.19 

1.70 

6.27 

± 

0.14 

3.40 

± 

0.06 

2.33 

± 

0.05 

14.01 

± 

0.17 

2.00 

4.67 

± 

0.14 

3.08 

± 

0.06 

2.36 

± 

0.06 

12.45 

± 

0.17 

2.40 

3.48 

± 

0.11 

2.54 

± 

0.06 

2.20 

± 

0.08 

10.59 

± 

0.15 

2.90 

2.22 

± 

0.10 

2.11 

± 

0.07 

1.93 

± 

0.09 

8.52 ± C 

.16 

3.40 

1.44 

± 

0.09 

1.58 

± 

0.07 

1.64 

± 

0.08 

6.72 ± C 

.15 

4.10 

0.95 

± 

0.08 

0.98 

± 

0.07 

1.29 

± 

0.06 

N/A 




TABLE 11. Fourth order truncated moments of the neutron 
F 2 structnre function based on the BONuS data for the three 
resonance regions studied, as well as for the whole region 1.3 < 
< 4 GeV^. For = 4.10 GeV^ the moment for the whole 
region was not computed due to the lack of data in the highest 
region. 



lOOO^Ms 

Q'^ (GeV") 

first 

second 

third 

whole 

1.00 

4.39 

± 

0.18 

0.60 

± 

0.01 

0.20 

± 

0.01 

5.28 

± 

0.18 

1.20 

4.19 

± 

0.10 

0.82 

± 

0.01 

0.30 

± 

0.01 

5.45 

± 

0.10 

1.40 

3.79 

± 

0.09 

0.98 

± 

0.02 

0.39 

± 

0.01 

5.38 

± 

0.09 

1.70 

3.24 

± 

0.08 

1.09 

± 

0.02 

0.52 

± 

0.01 

5.17 

± 

0.08 

2.00 

2.62 

± 

0.08 

1.13 

± 

0.02 

0.62 

± 

0.02 

4.82 

± 

0.09 

2.40 

2.12 

± 

0.07 

1.06 

± 

0.02 

0.68 

± 

0.02 

4.41 

± 

0.08 

2.90 

1.45 

± 

0.07 

1.00 

± 

0.03 

0.71 

± 

0.03 

3.77 

± 

0.08 

3.40 

0.99 

± 

0.07 

0.82 

± 

0.04 

0.67 

± 

0.03 

3.14 

± 

0.09 

4.10 

0.68 

± 

0.06 

0.56 

± 

0.04 

0.60 

± 

0.03 

N/A 



TABLE III. Sixth order truncated moments of the neutron 
F 2 structure function based on the BONuS data for the three 
resonance regions studied, as well as for the whole region 1.3 < 
< 4 GeV^. For = 4.10 GeV^ the moment for the whole 
region was not computed due to the lack of data in the highest 
region. 
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obtained from the BONuS data to the neutron moments 
calculated over the same x range based on the ABKM 
QCD fit [^. These truncated moments are shown as 
a function of in Fig. (closed circles) for M 2 . For 
the M 2 the earlier, model-dependent results from m 
are also shown (open circles). The four panels corre¬ 
spond to the four invariant mass regions previously de¬ 
fined. The computed ABKM moments include higher 
twist effects and target mass corrections, which allows 
for a direct comparison of the present results with the 
studies from Ref. m- The figure also presents the ratio 
of truncated F 2 neutron moments computed using the 
CTEQ-Jlab (CJ12) parameterization [21] with respect 
to the same ABKM moments. The explicit inclusion of 
higher twist terms in these parameterization allows for 
extending the region of validity into the large-x region. 
The interplay between higher twist and target mass cor¬ 
rections and their influence on the PDFs at large x is 
studied in depth in Ref. [3D]. In addition, the CJ12 pa¬ 
rameterization includes nuclear effects in part to improve 
the neutron, and hence d-quark PDF. 

To study the possible influence of target mass correc¬ 
tions and higher twists the CJ12 model was used to ob¬ 
tain both the leading twist (LT, thin dashed line) and 
the target mass and higher twist corrected (TMC-I-HT, 
thick solid line) predictions. The differences between the 
two models including higher twist and target mass correc¬ 
tions presented are within the uncertainty of the data. In 
the resonance region, target mass corrections and higher 
twist effects are expected to be sizable and so it may not 
be surprising that agreement of the data moments with 
leading twist parameterizations is improved by the inclu¬ 
sion of HT and TMC effects. However, the ABKM and 
CJ12 parameterizations extracted TMC and HT correc¬ 
tions from deep inelastic scattering data at higher Q^. 
Therefore confirmation of duality indicates that these 
contributions are similar in the resonance region and in 
the DIS. An alternative possibility is that the PDF hts 
have a very large uncertainty at high x and perhaps need 
to increase the d{x) values to better describe the data. 

In Fig. j^the first W range corresponding to the A res¬ 
onance displays the same behavior, but is 20-40% 
larger than the model expectation. This underscores 
the likelihood that the PDF-based fits underestimate the 
large x regime, to which this W region corresponds. 

Although the BONuS results reported here are in good 
agreement with the earlier data extracted from deuteron 
using modeling (open circles) [13], and both are consis¬ 
tent with unity in the second and third resonance region, 
the BONuS data are consistently higher in the first res¬ 
onance region. Although all local duality studies on the 
neutron are sensitive to binding effects, Fermi motion, 
and hnal state interactions, obtaining neutron structure 
functions by subtracting (smeared) proton from deuteron 
data is especially sensitive to the nuclear models em¬ 
ployed at the largest x (i.e. smallest W) values. There¬ 
fore, this discrepancy could be due to either a difference 
in the free neutron/proton A resonance excitation or a 



FIG. 4. The ratio between the M 2 truncated moments for 
the neutron structure function F 2 obtained from the BONuS 
data and the same quantity obtained from the ABKM pa¬ 
rameterization [ 27 ] (closed circles). The ABKM calculation 
includes both higher twists and target mass corrections. The 
open circles represent the model-dependent results obtained 
by [HJ. The lines represent comparisons of the ABKM with 
the CJ12 parameterizations |29| : thick solid line for the case 
when TMC and HT effects are included in both models, thin 
dashed line when the CJ model is evaluated only for leading 
twist. 


manifestation of nuclear theory uncertainties. Interest¬ 
ingly, the BONuS/CJ12 ratios assuming no higher twist, 
would fall consistently below unity, as evidenced by the 
dashed line. 

Moments in the resonance region display the same 
dependence as the PDF predictions, i.e. the ratio of the 
two shown in Fig. is independent of Q^. This is a re¬ 
markable confirmation of one aspect of duality, i.e. that 
the large x and lower region can be described by the 
same evolution equations - displaying the same dynamics 
- as the higher , larger W, scattering regime of single 
quark deep inelastic scattering. This is true even for the 
very largest x, lowest W, area. Here, the uncertainties 
on the PDFs are large and could easily account for the 
~ 20% strength difference, but not the Q^-dependence. 

Local quark-hadron duality in the proton F 2 struc¬ 
ture function has been studied extensively [IIS], show¬ 
ing possible violations mainly in the first resonance re¬ 
gion. To compare the proton with the neutron results, 
truncated F 2 moments were constructed for the proton 
using the ABKM model to calculate the deep inelastic 
moment and a global fit m for the resonance region. 
Fig. [5] shows the ratio of the neutron-to-proton trun¬ 
cated moments (closed circles) as a function of for 
the four resonance regions. The line shown is the ratio of 
neutron and proton truncated moments calculated using 
the ABKM (including HT and TMC) model. The model 
agrees well with the data in the second and third reso¬ 
nance regions but substantially underestimates the data 
in the first resonance region. This can be indicative of a 
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FIG. 5. The ratio of neutron to proton F 2 truncated moments, 
M 2 /M 2 for the four regions studied. The closed circles rep¬ 
resent the BONuS results, while the solid line corresponds to 
the ABKM QCD fit prediction including HT and TMCs. 


fundamental difference between the proton and the neu¬ 
tron at the very largest x values corresponding to the A 
resonance region. Alternatively, this can be the result of 
the fact that QCD fits for either proton or neutron are 
not well constrained at large a;, which is the kinematic 
range of the first resonance region. 

In conventional parameterizations of parton distribu¬ 
tion functions the d/u ratio is often assumed to go either 
to zero or infinity as x —?> 1, depending on the parame¬ 
terization of the d PDF [28]. At x above 0.8 the uncer¬ 
tainties on the PDFs, especially d, are considerable, due 
to deuteron nuclear model corrections and to the lack of 
large x data. The ratio d/u, and correspondingly F 2 /F 2 , 
at large x is however of fundamental interest [32l - l34| and 
there are numerous QCD-based theoretical predictions 
for this quantity. All require d{x) < u{x) and hence 
F 2 < F 2 for a given x and . This is in contrast to 
naive predictions for the A resonance. For the latter, 
given that the proton and neutron transitions to the A 


are isovector, the resonant contributions should be iden¬ 
tical, which would lead to an ratio close to unity. 

It is perhaps not surprising that the resonance behavior 
dominates in this W region because the A region contains 
only one resonance with less continuum background than 
the other regions. 

Similar tension between the resonance data and DIS/ 
pQCD-derived moments could exist in the second and 
third resonance regions as well. Assuming a dominance 
of magnetic couplings, as some quark models do [551137j . 
the proton resonance data should overestimate the DIS 
prediction in these regions due to the odd-parity reso¬ 
nances, such as the prominent spin—1/2 octet, resulting 
in a measured neutron-to-proton truncated moment ra¬ 
tio that should fall below the corresponding DIS curve. 
As shown in Fig. I^the results presented here do not sup¬ 
port this argument, and duality is apparently the pre¬ 
dominant effect for regions above the A. 


IV. CONCLUSION 

In conclusion, this paper investigates local quark- 
hadron duality in the neutron structure function based 
on data obtained by the BONuS experiment at Jefferson 
Lab, which used a novel experimental technique to tag 
the spectator proton in a deuterium target and thereby 
create an effective neutron target. This technique pro¬ 
vides smaller systematic uncertainties than earlier stud¬ 
ies that relied on the subtraction of smeared hydrogen 
data and nuclear modeling from deuterium. Truncated 
F 2 structure function moments were compared to PDF 
fits based on pQCD and largely deep inelastic scattering 
data, as well as to similar truncated moments obtained 
for the proton. The results indicate that quark-hadron 
duality holds dynamically everywhere, as well as in sum 
for the second and third resonance regions. This data 
provides further confirmation that duality, still not well 
understood, appears to be a fundamental aspect of nu¬ 
cleon structure. 
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